SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 

PHYSICS 

NANOSCIENCE AND 
TECHNOLOGY 

CONDENSED-MATTER PHYSICS 

THEORY AND COMPUTATION 

Received 
30 September 2013 

Accepted 
3 March 2014 

Published 
17 March 2014 



Correspondence and 
requests for materials 
should be addressed to 
U.S. (udo. 
schwingenschlogi® 
kaust.edu. sal 



Superior thermoelectric response in the 
3R phases of hydrated NaxRh02 



Y. Saeed, N. Singh & U. Schwingenschlogi 



Physical Science & Engineering division, KAUST, Thuwal 23955-6900, Kingdom of Saudi Arabia. 



Density functional theory is used to investigate the thermoelectric properties of the 3R phases of Na^RhOa 
for different Na vacancy configurations and concentrations. As compared to the analogous 2H phases, the 
modified stacking of the atomic layers in the 3R phases reduces the interlayer coupling. As a consequence, 
the 3R phases are found to be superior in the technologically relevant temperature range. The Rh d^^i-r^ 
orbitals still govern the valence band maxima and therefore determine the transport properties. A high 
figure of merit of 0.35 is achieved in hydrated Nao.83Rh02 at 580 K by water intercalation, which is 34% 
higher than in the non-hydrated phase. 

Thermoelectric materials are required for power generation from waste heat as well as for refrigeration. 
However, the practical application is limited by the low efficiency of thermoelectric devices as compared 
to traditional fuel power generators and compressor-based refrigerators'. A high performance thermoelec- 
tric material has to be a good electrical and poor thermal conductor and, at the same time, has to possess a high 
Seebeck coefficient^. The efficiency is determined by the dimensionless figure of merit zT = oS'TIk, where a is the 
electrical conductivity, S is the Seebeck coefficient, T is the temperature, and k is the thermal conductivity. 
According to the classical semiconductor theory, heavUy doped semiconductors with carrier concentrations of 
10" to 10^' cm"^ in general are promising candidates for high zT materials\ 

The discovery of a high Seebeck coefficient (S ~ 100 jiV/K at 300 K) in Na^Co02 has attracted considerable 
attention about a decade ago'''^. It has been predicted that the two-dimensional structure with alternatively stacked 
Na and C0O2 layers plays an important role for the high zT of Nao.5Co02''. Later, Kuroki et aU have argued that 
the peculiar pudding mold shape of the valence band of Naj.Co02 is responsible for the high zT (with high a), see 
also*. Interestingly, Takada et aU have demonstrated that the compounds can be readily hydrated to form 
Na^Co02 ■ 7H2O, maintaining the C0O2 layers, but with a considerably expanded crystal c-axis to accomodate 
the intercalated water. The electronic structures of hydrated and non-hydrated Na;cCo02 have been studied by 
first principles calculations in Ref. 10. 

Analogous compounds with Rh in place of Co are also good thermoelectric materials, since the octahedraUy 
coordinated ions Rh^* (4ii'') and Rh** (4d^) show similar electronic structures as Co'* (3d'') and Co"** (3(i^), but 
exhibit reduced correlation effects" Experiments could not verify any magnetic ordering in Na^Rh02'^. In 
general, the Rh 4d orbitals are spatially more dispersed than the Co 3d orbitals, which may result in an increased 
mobility of the charge carriers. Varela et al."' have synthesized Na^Rh02 as solid solutions composed of nano- 
metric particles of the a-NaFe02 structure type (space group R3m, no. 166) for x = 0.70 and 0.85. The 3R 
structure of Na;cRh02 contains three layers of RhOe octahedra per unit cell. Using the hexagonal setting of the 
R3m space group, the c-axis of the unit cell has a length of typically 18 A, while the length of the a-axis strongly 
depends on the specific transition metal cation. 

It is important to notice that the layered oxides LiRh02, NaRh02, and KRh02 can form hydrated, i.e., water 
intercalated, phases" with an elongated c-axis. For Na;rRh02 it has been shown that the water intercalation and 
the elongated unit cell result in a decrease of the resistivity'". Accordingly, in Ref 19 a colossal thermoelectric 
powerfactor has been demonstrated theoretically for K7/sRh02 in the case of hydration, which is nowadays 
triggering various experimental efforts^""". From a practical point of view the results reported on K7/8Rh02 are 
limited by the fact that the optimal thermoelectric performance is achieved below usual operation temperatures". 
In this context we propose to switch from the 2H phase (two RhOe layers per unit cell with a hexagonal structure 
of space group P63/mmc, no. 194, and c ~ 12 A) to the 3R phase, as one may speculate that the thermoelectric 
properties are modified by the different stacking (larger band gap in the 3R phases^^) . Accordingly, we study in the 
following the thermoelectric properties of both non-hydrated (c = 15.54 A) and hydrated (c = 20.85 A) 3R- 
Na;cRh02, for which we first determine the lowest energy structures for x = 0.75 and 0.83. Afterwards we analyze 
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the effects of hydration on the transport properties in the temper- 
ature range from 300 K to 800 K and demonstrate the great potential 
of the 3R-Na;cRh02 class of materials as thermoelectrics. The 3R 
phase turns out to be clearly superior to the 2H phase in the tech- 
nologically relevant temperature range. 

Results 

We start from the experimental lattice parameters a = 3.09 A and 
c = 15.54 A of NaRhOz'" and construct a 2 X 2 X 1 supercell. 
Afterwards we remove two Na atoms from different sites to obtain 
Nao.83Rh02. These configurations are named Al, A2, and A3. 
Similarly, we remove three Na atoms to obtain Nao.75Rh02 and name 
these configurations Bl, B2, and B3. All the structures under invest- 
igation are shown in Fig. 1. The formation energy is calculated for 
each configuration by the relation 

^Formation ^Defect ^Pristine ~t" 12(1 x) . 

Here Epristme is the total energy of a supercell of pristine NaRh02. In 
addition, fice/ecf ^nd are the total energies of the defective super- 
cell and of the isolated Na atom, respectively. The energy of an 
isolated atom is calculated by placing the atom in a cubic cell of 
20 A side length. 

The formation energies of configurations Al, A2, and A3 amount 
to 7.83 eV, 6.73 eV, and 6.65 eV, respectively, and those of config- 
urations Bl, B2, and B3 to 11.54 eV, 10.34, and 9.90 eV. These 
results show that for Nao.75Rh02 and Nao.s3Rh02, respectively, the 
A3 and B3 configuration has the lowest energy. The latter structures 
are characterized by maximal distances between the defects, which 
indicates that Na vacancies do not cluster in Na^Rh02, at least not in 
the two systems under investigation. Moreover, a minimum energy 
of 6.65 eV is required for creating two Na vacancies in NaRh02 and 
an energy of 9.90 eV for creating three vacancies. In the following, we 
will discuss only the lowest energy structures A3 and B3 unless 
otherwise mentioned. In pristine NaRh02 the Rh-O bond length 
amounts to 2.07 A and the thickness of the RhOe layer is 2.11 A, 



which is larger than the thickness of the CoOg layer (1.82 A) in 
NaCo02. Moreover, the bond angle /_ Rh— O — Rh = 96.5° is smaller 
than the bond angle ACo — O — Co = 98.35°. Next to vacancies the 
RhOe octahedra shrink slightly (the Rh-0 bond lengths decrease by 
0.02-0.05 A), because the interlayer coupling through the Na atoms 
is reduced and the RhOg octahedra thus become a bit more compact, 
which also influences the electronic and thermoelectric properties of 
the material. 

We determine the effects on the electronic band structure and 
corresponding density of states (DOS), for NaRh02, Nao.83Rh02, 
and Nao.75Rh02 in Fig. 2. NaRh02 has a indirect band gap of 
1.33 eV along the F-M high symmetry direction. The valence band 
mainly is formed by the Rh 4d tjg states and shows a total width of a 
bit less than 2 eV. This value is comparable to the bandwidth of 
NaCo02, even though Ad orbitals are more extended than 3d orbitals 
and thus usually exhibit higher bandwidths. The conduction band 
minimum belongs almost purely to the states. On the other hand, 
the valence band maximum, which is particularly important as it 
determines the consequences of hole doping, is due to the ci3z2_r2 
orbital. The electronic band structures obtained for configurations 
A3 and B3 are similar to the pristine case, except for the introduction 
of holes with densities of 8.8 • 10^" cm"' (A3) and 8.7 • 10'° cm"' 
(B3). Besides the characteristic peak approximately 0.2 eV below the 
Fermi energy, the DOS of prsitine NaRh02 shows dj,^i_y2 states 
down to —0.75 eV and between — 1 eV and —1.75 eV, which is 
similar for Nao.83Rh02 and Nao,75Rh02 except for the states at the 
Fermi level, which have different effective masses and consequently 
result in different Seebeck coefficients. 

The transport properties of hydrated and non-hydrated 
Nao.83Rh02 and Nao.75Rh02 are addressed in Fig. 3. We show the 
Seebeck coefficient and figure of merit for the temperature range 
from 300 K to 800 K. The Seebeck coefficient is positive in the whole 
range, indicating that the carriers are holes. We obtain room tem- 
perature values of S = 63 /(VK"' and S = 37 ^VK"' for non- 
hydrated Nao.83Rh02 and Nao.75Rh02, respectively. For hydrated 
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Figure 2 | Band structures and Rh 4d DOSs of NaRh02, Nao g3Rh02 (non-hydrated), and Nao.75Rh02 (non-hydrated). 



and non-hydrated Nao.83Rh02 the Seebeck coefficient behaves rather 
similarly as a function of the temperature, except for the high tem- 
perature range where hydration reduces S. Moreover, according to 
Figs. 3(c) and 3(d), zTat 300 K (room temperature) amounts to 0.16 
and 0.06 in the two systems. The former value is twice that found for 
Nao.85Co02 and Nao.88Co02 in Ref 5. Starting from almost the same 
value at 300 K, hydration enhances zT over the whole temperature 
range, such that at 800 K a value of about 0.4 is reached. At 580 K a 
maximal increment of 34% due to hydration is achieved. For 
Nao.75Rh02 the Seebeck coefficient increases at 300 K from 



37 /iVK"' to 51 ;(VK"' under hydration, where this increment of 
about 30% remains essentially constant over the whole temperature 
range. A similar effect is also found for zT. At 800 K again a value of 
about 0.4 is reached. 

The enhancement of S under hydration results in higher zT values. 
While the gross shape of the band structure changes little under 
hydration, see Fig. 4, the Rh 4d bands become flatter in a narrow 
region around the Fermi level. This fact causes the effective masses 
and therefore the resistivity to increase and results in a higher 
Seebeck coefficient. Moreover, the long c-axis of the hydrated phase 
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Figure 3 | Transport properties of non-hydrated and hydrated Nao.83Rh02 and Nao.75Rh02 as functions of the temperature, (a), (b) Seebeck coefficient 
and (c), (d) figure of merit. 



SCIENTIFIC REPORTS | 4 : 4390 | DOI: 1 0. 1 038/srep04390 



3 




causes the Rh 4d states to shift down in energy towards the Fermi 
level, because the Rh-O overlap is reduced (increased bond length of 
2.2 A). In order to study possible disorder of the vacancies we have 
also calculated the Seebeck coefficient (not shown) of the higher 
energy configurations Al, A2, Bl, and B2. For A2 and B2 the results 
are close to the values of the respective lowest energy configurations, 
because both have the vacancies in different RhOe layers. In contrast, 
concentration of vacancies in the same RhOg layer results in a sig- 
nificant reduction of S, which is, however, highly unlikely as certified 
by the previously listed formation energies. 

Discussion 

The Na;jRh02 oxides are found to form a new class of materials with 
exciting thermoelectric features, even outperforming the 2H phases 
of the Kj^Rh02 system. In the latter the optimal thermoelectric per- 
formance is achieved at low temperature, whereas the modified 
stacking of the atomic layers in the novel 3R phases of Na;tRh02 
results in a reduced interlayer coupling and, in turn, in a dramatically 
enhanced thermoelectric response in the technologically relevant 
high temperature range. We have also demonstrated that Na vacan- 
cies in Na;cRh02 avoid clustering and that the RhOe octahedra are 
modified depending on the amount of Na deficiency. Analysis of the 
induced changes in the DOS close to the Fermi level indicates that 
mainly Rh^'*''' d^^2_^i states control the transport properties of the 
compounds. A high zT value of 0.35 at 580 K is achieved in hydrated 
Nao.83Rh02 due to the enhanced effective mass of the charge carriers. 
In general, zT can be further increased by reduction of the Na 
vacancy concentration to increase the resistivity. 

Methods 

We study the effects of Na deficiency and hydration on the band structure of NaRh02 
using density function theory as implemented in the WIEN2k package"^"*. The 
exchange-correlation potential is parametrized in the generalized gradient approxi- 
mation^*. We fuUy optimize the structure of Na;cRh02 for x ~ 1, 0.83, and 0.75, 
considering 2X2X1 supercells, until the atomic forces reach values of less than 
10 mRyd/iJe- For the wave function expansion inside the atomic spheres a maximal 
angular momentum of ii^y^ax — 12 is employed and a plane- wave cutoff of R„,iK^a^ — 7 
with Gniax ~ 24 is used. Self- consistency is assumed when the total energy variation is 
less than 10"* Ry. The transport behavior is investigated using BoltzTraP'^^, which 
employs a methodology that previously has demonstrated to yield quantitatively 
accurate values for the thermoelectric properties of metals and doped semi- 
conductors^''"^^. It is based on semi-classical Boltzmann theory within the constant 
relaxation time approximation, i.e., the scattering rate is assumed to be independent 
of momentum and energy. This approximation allows the thermopower to be directly 
calculated from the band structure as a function of the carrier concentration and 



temperature""*. We note that it has been demonstrated for layered oxides isostructural 
to the present compound that at elevated temperatures the lattice part of the thermal 
conductivity can be neglected as compared to the electronic part"". We use a k-mesh 
of 7 X 7 X 2 points for the band structure and density of states (DOS) and a dense k- 
mesh of 18 X 18 X 18 points for the thermoelectric calculations. 
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